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Abstract 

A new one-dimensional fermion model depending on two independent in- 
teraction parameters is formulated and solved exactly by the Bethe ansatz 
method. The Hamiltonian of the model contains the Hubbard interaction and 
correlated hopping as well as pair hopping terms. The density-density and 
pair correlations are calculated which manifest superconducting properties in 
certain regimes of the phase diagram. 
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The study of strongly correlated fermion systems has attracted considerable at- 
tention during the last decade. In particular the discovery of high temperature su- 
perconductivity [1] has led to the investigation of various electronic mechanisms for 
an explanation of this phenomenon. In one space dimension several new integrable 
models have been found exhibiting different physical behaviour [2-9]. Notably the 
model with correlated hopping [5], the anisotropic tJ model [7], and models with 
pair hopping processes [8] show dominant superconducting correlations. 

Due to integrability many properties of these models can be found nonpertur- 
batively by an exact solution on the basis of a Bethe ansatz [10-13] and conformal 
field theory [14,15]. On the other hand integrability imposes severe conditions on the 
interaction parameters of a model. Usually, integrable cases are restricted to discrete 
points or depend on one parameter only. Therefore, the study of competition of 
different interaction terms is difficult. 

In this letter we report on a new integrable two parameter model of strongly 
correlated electrons unifying several previously known integrable systems such as the 
correlated hopping model and the tJ model. The phase diagram is very rich with 
a crossover from a regime with dominating density- density correlations to a regime 
with dominating (superconducting) pair correlations. 

Our starting point is a one-dimensional Hamiltonian with correlated hopping 
terms, with the Hubbard on site interaction, and with pair hopping processes 

r 1 1 

H= - I] (4aCj+ia + c++i^Cjv)exp --(r7-a7)nj-_^- -(77 + a7)nj+i,_^ 



+ E [Unj^riji + tp (4c+^cj+i|Cj+i| + h.c.)] , (1) 
j 

where j denotes the sites and we use standard notation for fermion operators. Tight 
binding Hamiltonians of this type have been considered since a long time [16,17] 
and have been studied as effective one-band Hamiltonians for the description of 
cuprate superconductors, see [18,19] and references therein. In (|I|) we have included 
an anisotropy 7 in the correlated hopping term and also a pair hopping term tp. 
Obviously, the integrable Hubbard chain [11] is a special case of (|l]) with 77 = 7 = 
tp = 0. For rj = ±7 and tp = U = we recover the correlated hopping model [5]. 
The exact solution for the particular case 17 = In 2, 7 = was studied in [9,20]. 
Here we shall establish the integrability of (|1|) under the condition 

tp = U/2 = e [2e-''(cosh r] - cosh 7)] e = ±1, (2) 

leaving us with two continuous parameters rj, {\r]\ > > for hermiticity) , and the 
discrete sign e where e = ±1 corresponds to repulsive and attractive pair processes, 
respectively. The exact solution for the eigenstates and eigenvalues of Hamiltonian 
(0) can be obtained from a Bethe ansatz. The structure of the Bethe ansatz equations 
follows from the solution of the two particle scattering problem. The nonvanishing 
elements of the 5* matrix are 



1 



^12 (A) 



^21 (A) 



sin A 



sin(A - 
zsinh7 



sin(A — i^) 
isinh7 

sm(A — 27) 



(3) 



where A = Ai — A2, and the \j are suitable particle rapidities related to the momenta 
of the electrons by 

^ ' I vr — 9(A, a), e = —1, 
with the function O and the parameter a defined by 



(4) 



6(A,a) = 2 arctan (coth a tan A) 
sinh + 7) 

sinh 5(?7 - 7) 
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(5) 



A necessary and sufficient condition for the applicability of the Bethe ansatz are the 
Yang-Baxter equations [10,21]. These equations are satisfied by (|^), i.e. by the S 
matrix of (|I|) under the condition (^). 

The Bethe ansatz equations are derived following the standard procedure by im- 
posing periodic boundary conditions. Each state of the Hamiltonian is specified by 
a set of particle rapidities A^, j = 1,..., N, and a set of spin rapidities Aq,, a = 1,..., 
N^, where is the total number of electrons and A^| the number of down spin elec- 
trons. All rapidities within a given set have to be different, corresponding to Fermi 
statistics. They have to satisfy the Bethe ansatz equations 



sin (Aj — ia) 
sin (Aj + ia) 
-iX- sin {Aq - Xj + i-f/2) 
j}^ sin (A„ - Xj - ij/2) 



N 

n 



I sin {Xj -Aa + i-f/2) 
\ sin {Xj — Aq, — ^7/2) 



-pi sin {Aq - A^ + 27) 
j-J-^ sin {Aa - Ap- i-f) ' 



(6) 



where L is the number of lattice sites. The total energy and momentum of the system 
are given in terms of the particle rapidities Xj as 



N 



E 



—2 ^ cos kj 

3= 

N 

2^E 



cosh 2a 



sinh 2a 



cosh 2a 



cos 2Aj J 



N 



E MA, 



(7) 



Next we discuss the physical properties of the model. First consider the limit 
7] — s> —00, e = +1. In this limit double occupations on the same site are excluded 
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and the remaining dynamics is identical to that of the supersymmetric tJ model 
[2,3,7] as can be seen from the Bethe ansatz equations and also from a canonical 
transformation of the Hamiltonian 



generated hy A = e l''l/^(c/ — c?+) where 



d= E {-l)'-'"'Pl+l,rr4+l.Cl.{l-ni^-.) 
i<m,(T=itl 



(9) 



Pa^b denotes the cyclic shift operator acting on the sites a, a + 1, 6, moving 
the state at site a to site h. In first order of A the effective Hamiltonian (H) does 
not involve interactions describing transitions between states with doubly occupied 
sites and those without doubly occupied sites. In the limit rj ^ —oo the resulting 
Hamiltonian is that of the supersymmetric tJ Hamiltonian. A more detailed account 
will be given in a separate publication. 

Secondly, a particle-hole transformation T together with a sublattice rotation 
(cjo- {—lycja) applied to H{ri) (|l|,|^) yields a resulting Hamiltonian of the same 
form, namely 

TH{t])T-^ = e'^Hi-T]) + (L - N)U. (10) 

Therefore we can restrict the further investigation to one specific sign of t]. We choose 
7] > corresponding to repulsive (physical) correlated hopping of the electrons and 
positive values for a (^. 

For the groundstate the magnetization is zero, Ni = N/2, and the second set of 
equations in (^ can be solved for the in terms of the Xj similar to the treatment in 
[5] (Appendix of the second paper). Inserting the result into the first set of equations 
in (^) and introducing the density function p(A) for the distribution of Xj in the 
thermodynamic limit, we obtain the linear integral equation 



27rp(A) 



2sinh2a /" ^ , ^ , 

— 7-^ ZT- + hp{X- n)p{^)d^, 

cosh 2a — cos2A Ji 



where 



2E 



-n7 



COS 2nA 



^^-^ cosh n7 
In order to minimize the groundstate energy 



^ = 2e 
L Ji 



the integration interval I in 



cosh 2a 



sinh 2a 



cosh 2a — cos 2A 



p{X)dX, 



(12) 



(13) 



|^|r^,|T^) has to be chosen symmetrically around (/ = 
[-K, +K]) or 7r/2 (/ = [7r/2 - fsT, 7r/2 + /T]), for £ = +1 or -1, respectively. The 
parameter K is determined by the subsidiary condition for the total density p = N/ L 
of the electrons 

'p{X)dX = p. (14) 
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Results for the groundstate energy per lattice site Eq obtained by numerical 
solutions of these equations are presented in Figs. 1 and 2 for various values of r] and 
7. For large values of the correlated hopping parameter 77 the processes described by 
Hamiltonian (|l]) simplify drastically. The only single particle process allowed is the 
hopping of electrons from singly occupied sites to empty sites with hopping rate t = 
— 1. In addition, pair hopping processes occur from doubly occupied sites to empty 
sites with parameters tp = U/2 = e. For e = +1 the groundstate consists of only 
single electrons for densities < p < 1, with dependence of Eo on p as for free spinless 
fermions {Eq = — ;^sin7rp). Upon further filling doubly occupied sites are created. 
However, the system freezes in as all hopping processes are suppressed. All states 
with arbitrary distributions of singly and doubly occupied sites are eigenstates of the 
Hamiltonian where the energy is just given by the Hubbard interaction. Therefore, 
Eq depends linearly on p for 1 < p < 2, see Fig. la). For the case e = —1 the 
formation of doubly occupied sites is favoured for all densities. The dynamics of the 
pairs is described by free spinless fermions {Eq = — ;^ sin7r| — p) for all densities p, 
see Fig. 2a). 

We next study the long distance behaviour of correlation functions by means of 
finite size studies and application of conformal field theory. For details of this method 
the reader is referred to [14,15,5]. There are two different kinds of excitations in model 
(0,11). One type of excitations corresponds to gapless particle-hole excitations, i.e. 
a redistribution of the Xj parameters, the second type consists of spin excitations 
with gap. The latter excitations become gapless for 7 = which situation will be 
studied in a separate publication. The central charge of model (|I],0) is c = 1. The 
long-distance behaviour of the density- density and pair correlations is given by 

(p(r)p(O)) ~ p^ + Air-^ + A2r~° cos(2fcj.r), (15) 



and 

(c+^c+^coiCot) ~ 5r"^, (16) 

where 

2kF = vrp. (17) 

The exponents a and (3 describing the algebraic decay are calculated from the dressed 
charge function ,^(A) which has to satisfy the integral equation 

27re(A) =27r + ^Vp(A-p)e(p)rfp. (18) 

We obtain 

« = 1//5 = e(i^)V2 (19) 

for £ = -1-1, and a similar equation with K replaced by 7t/2 + K for e = —1. 

The one-particle Green's function (c^^cqo-) shows exponential decay since hole 
excitations have a mass gap. Consequently the momentum distribution function 
{rikcr) is analytic in k, in particular there is no singularity at kp- 
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The model does not show finite off-diagonal long-range order. However, we ob- 
serve a longer range of the pair correlations in comparison to density-density corre- 
lations for certain regimes of interactions and particle density. Numerical results for 
the dependence of the exponent (3 on the density p are shown in Figs. 1 and 2. For all 
values of the interaction parameters rj and 7 there is a crossover from a regime with 
dominant density- density correlations (1 < /5 < 2) to a regime with dominant pair 
correlations {(3 < 1) at a "critical density" pc- For repulsive pair processes {e = +1) 
the size of regime "/3 < 1" is shrinking for increasing values of 1] and decreasing 
values of {pc ^ 2 for 7 — > 0), see Figs, la) and b). For attractive pair processes 
[e = —1) regime "/5 < 1" is growing for increasing values of 1] and decreasing values 
of 7 (pc for 7 —> and r] 00), see Figs. 2a) and b). 

In summary we have presented a new integrable model for strongly correlated 
electron systems depending on two free parameters. We have obtained exact results 
for the groundstate energy of the system as well as the critical exponents for cor- 
relation functions showing generically a crossover between regimes with dominant 
density-density correlations and pair correlations, respectively. 

After completion of this work we learnt about the formal proof of integrability of 
model d^) in the particular case 7 = [22,23]. 
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Figure captions 

Figure 1: a) Depiction of the density dependence of the groundstate energy per site Eq 
and the exponent (3 of the pair correlation function for values e = +1, 7 = 0.5, 
and r] = 0.6, 1, 2, 4. b) Similar to a) for e = +1, f] = I, and 7 = 0.1, 0.3, 0.5, 
0.7, 0.9. 

Figure 2: Similar to la) and b) but for e = —1. 
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